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ABSTRACT 

An investigation of n -p elastic scattering and single -pion pro- 
duction at 0.939 Bev/c is reported. Five-hundred and forty two-prong 
interactions in the 10-inch UCRL liquid-hydrogen bubble chamber have 
been classified. Emphasis is given the methods used in the kinematic 
classification of events. The total cross section for m -p interactions 
is found to be 53.3 + 2.4 mb. The elastic-scattering cross section is 
21.0 ۶ 1.1 mb. The elastic differential cross section displays a 
narrow diffraction peak in the forward direction. In the center -of- 
mass system, 35.6% of the elastic events have pion scatters into the 
backward hemisphere. 

In single-pion production the ratio of the number of r to n 
productions is 2.49. The momentum distributions of the charged 


particles do not correspond to the phase-space distributions predicted. 


lv 
















awaqpa ۲۳۸۱ نبا‎ 7 ۷ u ag“ Mua r 


ቺ መ < 
we, dnd 
2 Ú لام‎ 














E C የደ 









= 

CU o W 

w سا اس‎ smie = lı age 

፡‏ 11 ین ډب سسوم مر 0-0 اسه وص 

— — € مره قو د وسعسه مس‎ ዘ) ሐ. 
۳ ۱ 1 da ge dm n | 

—— — چب‎ 
en 1 
ييا‎ በመጨ | 






















E 


m -p ELASTIC SCATTERING AND SINGLE PION PRODUCTION 
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I INTRODUCTION 
An experimental study of elastic scattering and single-pion- 
production interactions of negative pions of 0.939 Bev/c momentum in 
hydrogen is reported here. 
Pion-proton interactions in this Bev region leading to associated 


Fae 


production have been investigated recently. Inelastic pion reactions 


not involving ''strange particles'' have been neglected except for early, 


4 


low-statistics experiments using diffusion cloud chambers and 
nuclear emulsions. : Elastic pion-proton scattering in this energy 
region, on the other hand, has received additional treatment using a 
propane cbe and, more recently, liquid hydrogen bubble 0 
The study described here is based on data obtained as a by- 
product of a recent associated-production experiment with the 10-inch 
hydrogen bubble chamber at the University of California Radiation 
Laboratory. ? Approximately 600 two-prong pion events have been 
measured at each of three (nominal) beam momenta: 0.927, 1.12, and 
1.24 Bev/c. The preliminary kinematic analysis of the 0.927 -Bev/c 
data is reported here. Analysis of the 1.12- and 1.24-Bev/c momenta 
data will be reported later. 
The following three reactions are considered: 
T + م‎ =+ DES 
ፐ ሰ زک‎ 
ፐር ቄር ፒክ 
The decision to limit the kinematic analysis to the three listed reactions 
assumes that the frequency of events with the production of two or more 


secondary pions 1s negligible. The rarity of four-prong events, 
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T +p—>T + و‎ + Y + 1 , supports this assumption. In principle, 
measurement of the charged-particle tracks in a magnetic field, and 
conservation of momentum and energy, determine unambiguously 
whether one or more neutral particles are involved. In practice the 
identification of missing neutral particles is made more difficult by 
measurement inaccuracies. Classification of multiple -pion-production 
events requires well-established criteria as to what constitutes a good 
kinematic fit in single-pion production. An approximate upper bound 
for multiple pion production will be established later. 

This report treats the interpretation of the experimental 
results in only a cursory way. Specification of a single-pion-pro- 
duction event requires five intrinsic parameters, in contrast with the 
single center-of-mass scattering angle required for an elastic event. 
Consequently a detailed interpretive treatment of the data is relegated 
to the future. To keep the experiment in perspective, however, it is 
useful to list objectives of a study of elastic and single-pion-production 
events in this energy region. 

a. Elastic t -p scattering 
l. Determination of the size and opacity of the proton. 
2. Evaluation of the relative weights that the various angular - 
momentum states contribute to the interaction. 
Test of the validity of the dispersion relations (15% effect). 
4. Confirmation of the previously asserted backward-scattering 
peak (at higher energies). | 
b. Inelastic single-secondary-pion production 
l. Isotopic spin dependence of the interaction; the ratio of 

T ton production (a priori of equal phase space). 

2. Evidence for the existence ofa T= J= 1 nucleon isobar. ር 
3. Correlation of the inelastic pion production with associated 


production. 
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II. EXPERIMENTAL PROCEDURE 

The r -p interactions analyzed in this experiment were selected 
from a group of photographs taken in the 10-inch-diameter hydrogen 
bubble chamber. The chamber was provided with a 10,900-gauss 
magnetic field. Details of the hydrogen bubble chamber used have 
been described. i The exposure was the codperative effort of the 
Alvarez group, and was primarily for the purpose of nvestigating 
associated production of strange particles at 0.927 Bev/c nominal 
beam momentum. The t -beam trajectory is discussed elsewhere. E 
About 24,000 Bevatron pulses were photographed during the run. 
Approximately 4,800 of these were used in this analysis. All film 


was scanned by the professional scanning staff. 


111. SCANNING EFFICIENCY AND ACCEPTANCE CRITERIA 
A. Scanning Efficiency 

l. General 

At each transit of the m beam through the bubble chamber, two 
stereophotographs of the chamber were taken. After processing, the 
35-mm film was scanned for all events of interest. Pion interactions 
were recorded, according to frame number and type, by means of 
check marks in appropriate columns of the scan sheets. Events 
entered on these scan sheets were subsequently examined in detail, 
sketched, and measured. The 13 rolls of film used in this experiment 
were scanned twice. Several factors which reduce the probability of 
seeing an event are discussed, together with their compensating 
corrections. 
2. Depth Dependence 

An event with its vertex near the top or bottom glass plates of 
the chamber is likely to be missed because one or both of the emergent 
prongs can be very short. An analysis of this chamber depth dependence 
was based on a comparison of the vertical distribution of 575 event 
vertices and the vertical distribution of 1126 beam tracks. The chamber 


was considered as divided into five horizontal regions. The median 
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plane of the central region, C, coincides with the median plane of the 
bubble chamber; depth of region C is 3.2 in. Regions A and E are the 
extreme regions, extending to the top and bottom glasses, respectively. 
Regions B and D lie between the central region and regions A and E, 
respectively. These divisions are each 0.8 in. in depth. 

The depths of beam tracks in every 50th frame were measured. 
The vertical position of each event vertex was obtained by measure- 
ment on the precision measuring engine. The ratio of event density 
to track density within each horizontal section of the chamber is a 
measure of scanning efficiency within that region. These ratios are 
as follows. 

Table I 


Chamber Depth Dependence 


Region Events(N) Beam Tracks(M) Ratio (N/M) 
A 2 2 100 7 
B 3 12 0:25 3 
ር 563 1082 95 9 
D ር 22 9 
E l 8 papa 

Total 575 1126 0.51 


Obviously scanning efficiency for events in depth region E is poor. The 
ratios for allother regions agree within statistics, and only the events 
found in these four regions were considered. No correction to scanning 
efficiency as a function of chamber depth was made. 
3. x-y Dependence 

A reaction that occurs near the chamber wall may also be 
difficult to detect because of short tracks. Aside from this consideration, 
momentum and angle measurements become less accurate as the length 


of track available for measurement decreases. Hence, long tracks 
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were desired. A central area was superimposed on one stereo view at 
the stereo projector, using chamber fiducials for positioning; only 
those events whose vertices projected into the area were considered 
for measurement. This criterion not only increased the scanning 
efficiency, and the probability of having long tracks, but it also made 
identification of events from the decay, interaction, or stopping of one 
of the emergent prongs more likely. 

In order to detect any variation in scanning efficiency within 
the acceptance area, this area was divided into five equal regions, as 
shown in Fig. 1. The distribution of 600 projected vertex positions 
was compared with the distribution of 8562 beam tracks crossing the 
area. Ali tracks in every tenth frame of the experiment were counted. 
On the histogram of beam track distribution, Fig. 2, the number of 
events within each subdivision of the acceptance area has been indicated. 
The ratio of event density to track density within each region of the 
acceptance area is a measure of the scanning efficiency for events 


within that region. These ratios are as foliows. 


Table II 


x-y Dependence 


Region Events (N) Beam Tracks(M) Ratio(N/M) 
1 91 1415 0.064 + .009 
2 122 1784 0.068 + .008 
3 2 21615 02072 2,007 
4 141 1871 0.075 + .009 
5 94 1377 0.068 + .009 
Total 600 8562 0.070 


These ratios indicate that scanning efficiency is uniform, within 
statistical accuracy, throughout the acceptance area, and that no 
dependence of scanning efficiency on position in this area need be 


considered. 
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Fig. 1. Projected view of bubble chamber, showing acceptance 
area. 
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Fig. 2. Projected x-y distributions within acceptance region. 
Histogram of distribution of 8562 beam tracks. (600 events) 





4. Orientation of Outgoing Prongs 

Events that have one or both prongs near the vertical are more 
likely to be overlooked than events with prongs more nearly horizontal. 
To study the dependence of scanning efficiency on orientation of the 
plane defined by the two prongs, a plot of the number of events as a 
function of the angle $ was constructed. Here $ is the angle between 
the horizontal plane and the plane containing the two prongs. This 
plot is shown in Fig. 3. The scattering should be isotropic in $, and 
the anisotropy observed is a measure of the dependence of scanning 
efficiency on event orientation. A correction factor of 1.18 + .03 is 
derived from the data illustrated in Fig. 3. The total number of 
events seen must be multiplied by this correction factor to account 
for the scanning-efficiency dependence on 6. 
5. Repeated Scans 

All film used in this experiment was scanned twice by two in- 
dependent observers. Within the acceptance volume there were 922 
events seen by both observers; 33 events were seen by the first scanner 
but not seen by the second; 27 events were seen by the second scanner 
but not by the first. The indicated efficiency of each scanner is 
97% + 1%. The total number of observable two-prong events, based 
on these scans, is 983 +1. The total number of two-prong events, 
after correction for scanning-efficiency dependence on $, 1s 1160 + 31. 

B. Acceptance Criteria 

Each two-prong event within the acceptance area was examined 
at a scanning table before being measured. Acceptance criteria, 
specified to reduce analysis time and measurement errors, and to 
fully utilize the automatic tracking capability of the measuring engine, 
were imposed on each event. These criteria were: (1) The vertex 
was distinct, 1. e. , no track not associated with the event fell within 
one track width of the vertex in either projected view. (2) The 
incident track direction was within 3° of the mean incident beam 
direction. (3) The incident track curvature was approximately the 
same as the mean beam curvature. (4) Not more than thirty beam 


tracks were visible in either view. 
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Fig. 3. Orientation of plane of emergent prongs. Folded 
distribution of 550 events. § = 0° = flat in chamber. 
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IV. TRACK ANALYSIS 
A. Spatial Reconstruction of a Track 

Each of the two views of an event accepted for measurement 
was projected on the screen of the precision measuring engine 
(Franckenstein). This servo auto-tracking projector takes a series 
of points along each prong of the event in each of the two stereo views. 
The solution of the optical problem and the determination of a set of 
points in chamber space from the Franckenstein measurements have 
been previously described. 

When the space coordinates X; Yo Z; (1-2 1,2,--- wor each 
of the n measured points along a track have been determined, a least- 
squares fit of a parabola is made to the projection of these points on 
the x-y plane.(Fig. 4.) The least-squares fit is made on the x.” and 
ሃ;- points following a translátion and rotation to the x'', y'! coordinates 
used to simplify the scaling problem for the digital computer. The 
x", y' coordinates are fitted to the parabola 

Yi. 2 2د‎ ® x, FC e 
(The event reconstruction described is programmed for the IBM-650 
computer.) 

Although the projection of the track on the x-y plane is a 
circular arc, a parabolic fit is made rather than a fit to a higher- 
order curve, because at the momenta measured the increase in 
accuracy 1s not sufficient to justify more elaborate programming. 

When the parabola has been determined, the root-mean-square 


deviation Yan of the measured points from the fitted curve is ex- 


8 2۱ 2 


xy n-3 ' 


pressed by 


G 





where n is the number of points measured, and 5. is the deviation 
of each projected point from the fitted curve. For n cox 7 is not 
meaningful, since a parabola can be fitted exactly to three points, and 
therefore no calculation of e is made. 

The azimuthal angle, between the tangent to the x-y projection 
of the track and the x axis, is given by 9 = 6, + oue The value 


R 
of በ” is obtained from the coefficients  b' and c! and the half-chord 
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Fig. 4. Space coordinate system and track projections into 
planes. 
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length R: 
R = x" - x. 

For each track there are two azimuthal angles, $; and ቀ፥› correspond- 
ing respectively to the beginning and end of each track. The value $. 
is used only with the incident beam track, for on this track alone is 
measurement directed toward the event vertex. 

To specify the length of the track, £, and the dip angle ^, the 
coórdinates are fitted to a straight line in the x''-z' plane. Although 


it is true that each track is a helix, and hence its projection on the 


x'-z' plane is sinusoidal, the momenta were sufficiently high to allow 


a least-squares fit to be made to the straight line و2‎ = d + ex". 
The root-mean-square deviation, defined by 
22م‎ 
2,0. V 
0 = 1 
Z n-2 ۱ 


is then obtained from the coefficients of the straight-line expression. 
As above, Š, and n are the deviations of the measured points (pro- 
jected into the x''z'' plane) from the fitted curve and the number of 
points measured, respectively. 

With the z dependence of the track now specified, the dip angle 
A and the track length £ are given by the expressions 
23 | 


tan À = - EE FR 


e. 
R ? 
Fitted end points are obtained from the least-squares solution. 

The final quantity of interest--the projected curvature of the 
track K--is determined from the fit in the x-y plane, according to 
K = 1/.3Hp, where p = Roe The magnetic field strength, H, is 
evaluated at the midpoint of the track. The momentum can be calculated 


from the curvature K and dip angle X by the expression p - (K cos (00 
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B. Measurement Errors 

The automatic tracking feature of the precision measurement 
projector makes positioning of the stage at each point operator -inde- 
pendent. This feature makes it feasible to accurately determine the 
errors associated with measurements of angle and curvature. 

Semi-empirical formulae for the standard deviation errors in 
the curvature K of a track, the azimuthal angle ቀዌ, and the tangent of the 
dip angle, tan ^, were constructed from statistical studies of meas- 
urements made both on beam tracks and on tracks of actual events. 

Errors fall into two categories: (a) internal-measurement 
errors, which reflect the accuracy of the measuring device and the 
intrinsic image quality of track; (b) turbulence errors, resulting from 
distortion of tracks by the turbulent liquid hydrogen during the time 
interval between occurrence of the event and the light flash. 
l. Internal Measurement Errors 

During measurement, track points are taken at approximately 
l -centimeter intervals along the projection of the track in the x-y 
plane. A fit to these measured points ES ے٦ ۰۶ء۰ ,1 ,90ز‎ 6+ 
equal standard deviations for the y setting for all points, is specified 


as follows: the ٠ quantity, 
n m à 
"”->፦ 3 (y- 2 ጻጓ] 
© [20 ال‎ 70 


is required to assume a minimum value for the "best" trial value of 
the parameters a; (i -0,1,:--m). The intrinsic accuracy of the 0: 
once determined, leads directly to the standard deviations Ou SE 
and UE , attributive to measurement technique. 

Let us consider the error in the dip tangent. If the z dependence 
of a track is stated as z - 50 + a) x, with x the distance along the 


projected chord length of the track, the dip tangent is given by 
fanen 2 02 dx = aj. The standard deviation is then 


NRG ° 
u 5 > ere" een رہ‎ በ 
Pp eee مس‎ "e 
inas: MT مت‎ 


هده مجو رون وسكت 





-14 


0 = G 


۲ 
tan À "61 
or 
0 = E zx G À (DSL) 
tan À R NN N 1 


where R is one-half the projected chord length, N is the number of 
points measured, and C,, is a correction factor for dependence on the 
number of points taken. 

The derivation of error in azimuthal angle follows a similar 
pattern. When a fit has been obtained in the x-y plane of the chamber, 
following coórdinate rotation, of the form 

y= + aX + a,x”, 


then examination of Fig. 4 shows: 


tan(,-,) = dy/dx| وی‎ 2 ۰ 
For a 10-cm track with momentum of 150 Mev/c, $: -OR = 14.7%, and 


this angle decreases rapidly with increasing momentum. Setting 


ር05(ቀ, -ቅኳ) = | allows the error in ቀ, to be len 


~ 4N3 cÊ 


0 = — ሀ : DU 
Q RN xy N | | 
On the basis of the equation for the fitted curve, as given above, 


the curvature of a track can be defined as 


K - (a^y/ax^) [1 * (ay/dx)^] E 


Evaluation of this expression at x = R gives K= 2 0 « where, by 
choice of coordinate axes, the PE contribution can be neglected. 
Propagating the error of fit in a, one can approximate the curvature 


error by 


Uu ل3‎ . (IV -3) 


K RNN xy N 
Both 6 and 07 are calculated for the fit achieved. From 


repeated measurements of individual events and measurements of 





10. 


beam tracks, mean values of 0.0041 cm and 0.0558 cm, respectively, 
2 with Eqs. (IV-1) to (IV -3), 
B = 305 cm/Bev/c where B = 1/(pK), based on the mean 10.9-kilogauss 


were found. Use of these o. and © 


magnetic field, results in the following approximate expressions for 


internal measurement errors: 


= 05 ۱ a 

o A A e (IVA) 
tan A RNN 
چ‎ ያ 77+70 (IV-5) 


ቀ RNN 


17. 


i 8.30 ይም ክ 
፲ኗ RNN N 


0 (IV -6) 
The ቅቆ dependence of لئ‎ reflects the use of both stereographic 
views in the determination of tan\\. _ 
2. Turbulence Errors 

The superheated condition required for track formation is 
obtained by rapid expansion of the liquid hydrogen prior to entry of 
each beam pulse. Because of turbulence following each expansion, 
tracks in the chamber are distorted, and curvatures obtained from 
photographs of these tracks can differ considerably from the curvatures 
of the actual ionizing particles. 

The momentum-error data of ہے‎ were taken as the basis 
for including turbulence effects within the total curvature error, Ok: 
Measurements of 386 beam tracks at 1.12 Bev/c momenta were made. 


The square of the standard deviation of curvature can be expressed as 


N 2 
l 1 1 
9. = N] » | AS ورد‎ | 9 (IV 7) 
1 Pi P / mean 


where es is the mean value of measured curvatures. Data were 
obtained by Cresti under the following conditions: (a) photographs were 
taken with 3 msec. light delay; (b) tracks were in the central region 

of the chamber (the exact form of turbulence is dependent on the 


position of the track with respect to the chamber expansion line); 
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(c) the average interval between measured points was 1.0 cm. Meas- 
urements were made on beam-track segments of varying lengths; three 


average track lengths were used to determine the R dependence of 


Du. The average lengths of the tracks considered were 6.1, 9.8, and 


19.8 cm. The three points of Fig. 5 were determined by use of Eq. (IV -7) 
and the Cresti measurements. 


A fit of these points is possible with a quadratic expression 
Ds. A T Bx + Guo where x - RU and R, the half chord length of the 


projected arc. Such a fit is indicated by inspection of the points. The 
correspondence between the functional form of the internal error in 
! 
curvature, O, in Eg (IV), and the last term of this expression made 
Hed 


possible a natural division of the total error, i.e., O e (G ) da ORE 
፻ ! K K K 


where Or and Ok are the turbulence and internal-measurement errors 
respectively. The dependence of C' on N, the number of points meas- 
ured on a track, was taken from the internal-error expression as 

C (N) = a (1 - x For flat tracks with points taken at 1.0-cm intervals, 
one has N =2R, hence the expression 


"ር - A* BR^ * 0.5 C! RÓ(R^ (1-0.75 R^ 


was fitted to the points indicated. The coefficients A, B, and C' 


1)] (1-8) 


required for this fit were 


A = 9.32 X j B = 0.174, 0 ۷٣٣ 


Repeated measurements of single events were also used to distinguish 
between internal errors and turbulence errors. 

In the kinematic analysis of events, Section V, the quantity 
. » where 5, = (x; - zu O., (i = 1,2,::: 9), is a measure of the 
shift made in the ith variable to obtain a fit to a designated 
interpretation. The distributions of Eb 2 and 55/0 ع‎ over a 
number of different events are a measure of the fidelity of the total- 
curvature error formulae when applied to the outgoing negative and 
positive track respectively. Histograms of these distributions are 
shown in Fig. 6. 

The distributions of 5:/0 e. over many different events 
complement the repeated-measurement technique, and these distrib- 


utions were used in the adjustment of the error-formulae coefficients. 
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Fig. 5. Curvature error from study of 386 beam tracks. 
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with stopping proton). 








aoe 


The error formulae as programmed for normal nonstopping tracks 
pping 





are | 
Ok - [A+B/R- Í RAE Ber on (IV -9) 
where 
A = 18.64 X 1077. 
B - 0.226, 
) 2027 ۱٢ IS not), 
R = + projected chord length (cm), 
N = number of points measured in View زا‎ 
o. =[0' 4 +0" ° 13 (degrees) (IV-10) 
$ $ ቀ 
where 
, 2 =] => 
D 0ت‎ ۲ N )1 ۰ 
UE - (0.188 R )* (A * B/R^); 
and 
Ou xs Sud. (1+2] sino 1)7?. (IV -11) 
"^ 2 


The o term in Eq. (IV-10) reflects the effect of turbulence 


t? 


b’ 


on azimuthal angle q. It can be shown that this turbulence error, O 
tt 
and the turbulence error in curvature, 0 x: are related by the 


expression 
! 


o (rad. | = (R/305)0 y. 
where 305 cm/(Bev/c) relates curvature and momentum in the magnetic 
field, and 0 y = (A +B /R2)2, as specified in Eq. (IV-9). 

For protons that stopped in the chamber, the momentum (and 
hence the curvature) was determined from the range -momentum 
formula, p = 0.141 r where £ is track length (cm). This | 
expression, based on the range-momentum curves for liquid hydrogen, 7 
defines the proton momentum much more accurately than is possible 


when curvature measurements are based on parabolic fittings. 
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Ihe range-momentum expression can be approximated by 
p “O07141 4 for determination of curvature error associated with 
measurement of a stopping proton. With half projected chord length, 
R, and dip-angle tangent, tan A, as independent variables, differen- 


tiation of this expression gives 


[fo 0 E 
1 R | zer? 
~ = — E ) 
DER cos] no TU | (Bev/c), 
where 4 = 2፪/ር05 እ. For "flat" tracks this expression reduces to 


Ok Ko, /4R. A value of 0.055 cm was taken as the error in one-haif 
chord length, Tp: This value was the standard deviation in range for 
444 TM ranges measured by Clark and Diehl, 7 usıng the same meas- 


uring engine. The approximation programmed was 


0 0.014 K/R. በ (IV 12) 


de 
Since the error in K for stopping protons is essentially 
internal-measurement error, 36 remeasurements of event No. 495,621 
were made to test this approximation. The arithmetic mean of the 
36 O'S for the stopping proton obtained by Ex. (IV-12) was 0.072 je 
The K, distribution for these remeasurements, shown in Fig. 7, has 
a standard deviation of 0.128 Vc The disagreement between 
the calculated error and the error associated with a typical event 
indicates the dependence of range-measurement accuracy on track 
characteristics. The measurements by Clark and Diehl were re- 
stricted to flat tracks with angles of intersection--at the vertex-- 
between 30° and 150°, and to the less heavily ionizing muonand electron 
tracks--that is, conditions not normally encountered with actual re- 
coil protons. These considerations suggested that the expression for 


error in K be broadened to 


gh s (ON. Bee (Iv 13) 


The typical accuracies of the two methods of measuring 
curvatures (momenta) are illustrated by event No. 494,808: 

track length (£) = 3.034 cm. ; 

\ =.58.59°; 


momentum by curve fitting, Po 7 121 د‎ 116 es 
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Fig. 7. Curvature of stopping proton; 36 repeated measurements 


of event No. 495,621. K, = 9.68/Bev/c; OK. = 0.128/Bev/c. 
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momentum from range, Pp, = 192.2 Mev/c; 

momentum from least-squares adjustment, Pr = 188.5 Mev/c; 

momentum-measurement error from O Eq. (IV-13), = =5.0 Mev/c. 
Here p. and P. differ from Pe the "best" value of the momentum, by 
2% and 36%, respectively. The adjustment of p, to obtain a "fit" was 
0.740 . The inaccuracy of p. reflects the short projected length of 
track (1.63 cm) available for the parabolic fitting. 
3. Incident-Beam Momentum 

The nominal beam momentum of 0.927 Bev/c was determined 
from a study of pion trajectories through the collimation system. 
The absolute value of the mean pion-beam momentum was determined 
from an analysis of 77 elastic events in which the proton stopped in the 
chamber. The nominal beam momentum, with a standard deviation of 
0.278 Bev/c, was used in the kinematic fit program discussed in 
Section V. This large, unrealistic error allowed unrestricted shifting 
of the beam curvature, Ki. in obtaining a fit to the intrepretation. In 
effect, no exact beam momentum was assumed. The adjustment 


parameter M 


where x, and x are respectively ''best fit'' and nominal values of 
Kı» ı 22 -0.278 Bev/c, indicates the adjustment of K} 
required for each event to fit the interpretation. The 51 histogram 


and 0 


for the 77 events is given in Fig. 8. The mean §, of -0.0414 + .005 
corresponds to a mean beam momentum adjustment to 0.939 € .006 Bev/c. 
The mean beam momentum, from a study of three associated production 
events of the type مہ‎ +p» A+ Ee and t + p> va obtained during 
this exposure, gives p, - 0.938 + .002 Bev/c. = A beam momentum 

Bev/c was used. The indicated beam spread,‏ 007. € 0.939 = را 

+ 0.007 Bev/c, is supported by a detailed study of associated pro- 
duction at 1.12 Bev/c. The beam trajectories through the magnetic 
steering magnets and through the collimation system were in all 
essentials identical with the 0.927-Bev/c experimental setup. The 
study of stopping protons, on the other hand, indicated only that the 
beam spread was less than + 0.035 Bev/c, a width that reflects the 


large intrinsic measurement error of the method. 
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Fig. 8. Adjustment of beam curvature for least-squares fit. 
Mean 3 = -0.414 + .005; 7 = 0 1727; 
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V. KINEMATIC ANALYSIS 
A. General 

In the past, 1t was possible for the physicists using cloud 
chambers to keep up with the data accumulation by "hand analysis " 
techniques. The bubble chamber, with its relatively high density, has 
greatly increased the number of interesting events observed, and 
therefore has stimulated the development of rapid analysis techniques. 
In this experiment an event was observed with approximately one out 
of every two Bevatron pulses, or 300 events were photographed each 
hour. The analysis of the events photographed during one 24-hour 
period of operation would require several man-years of "hand 


analysis. " 


Rapid precision analysis techniques are obviously required 
if the rate of analysis is to parallel the accumulation of data. 

All computations described in this report were carried out by 
use of the IBM type 650 digital computer. The procedure followed 
in the analysis of two-prong events is outlined in Fig. 9. The meas- 
urement and event-reconstruction phases of analysis have been dis- 
cussed in Section ÍV. 

B. Computation of Measurement Errors and Coplanarity 

The origin of measurement errors and the derivation of 
requisite formulae are discussed in Section IV. The computer is 
programmed to calculate the errors associated with all angle and 
curvature measurements. A section of this program computes the 
included angle between any two tracks; individual track errors are 
propagated to these included angles. 

The coplanarity of an event may be measured by the value of 
R = 81: (n, x راو‎ where ny: 82, ñ, are the unit vectors in the 
direction of the incident tm and the two outgoing prongs, respectively, 
at the vertex. Perfect coplanarity is, of course, indicated by R = 0; 
however, because of measurement errors, events that are actualiy 
coplanar may have values of R ደሽ መውን) coplanarity function F and its 
555 0-614۲06060707, Op, were used in establishing the degree of co- 
planarity of each event. Here R and F are related by 


R = F cos NJ cos A cos Az 
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Fig. 9. Flow diagram for analysis of two-prong event; computer 


time indicated. 








G 


where 
ps tanh, sin (6,-$,)+tan A sin(d,-$,) ttan Az sin($, -6,). 


Since tracks dipping or rising vertically in the chamber are not seen, 
cos À is never zero, and R = 0 if and only if F = 0. The error, Op: 
is determined directly from the measurement errors computed for 
each track. The quantity 7 0/60 indicates coplanarity within meas- 
urement errors. Figure 10 is a histogram of "27ا‎ |/ values for 
elastic and inelastic events. For the inelastic events a broad 
distribution of IF1/ 0, is expected, since the lack of coplanarity is 
determined by the direction and momentum of the neutral particles. 
C. Graphical Analysis of Inelastic Events 
1. Ellipse Plots 
In the following discussion, inelastic events 

N O, 
and 

T+p — a ar T tn 
are represented symbolically as 

(M,) + (M,) - (M,) + (M4) + (Mg). 


Subscripts on kinematic variables indicate the particle to which the 
variable applies. Unstarred quantities and starred quantities pertain 
to the laboratory and center-of-mass (c.m.) systems respectively. 

In the laboratory system that component of the momentum of the 
neutral particle, M,, that is perpendicular to the plane defined by the 
two charged particles introduces a deviation from coplanarity and 
prevents the use of a geometrically simple dynamical representation 
in this reference frame. Inthe c.m. system the three final-state 
particles are coplanar, and, in this frame, ellipse plots may be used 
in the kinematic analysis. The ellipse plots for analysis of the.two 
inelastic reactions are shown in Fig. 11. 

The properties of each ellipse are established by the Lorentz 
transformation from the barsy system, (primed quantities) in which 
the virtual mass of particles Mz and Mg is at rest, to the c.m. system 


(see Fig. 12). 
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Fig. 10. Coplanarity of 250 elastic events (solid line) and 250 
inelastic events (dotted line). 
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12. Ellipse parameters, with diagram of transformation 


Fig. 
from barsy to center-of-mass system. 
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Thus for particle No. 3, (p', = -p 
can write 
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p 3 COS De 
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a 5 . | 
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only the longitudinal component of the momentum -energy vector being 
transformed. The quantities y and n refer to the motion of the barsy 
relative to the c.m. system. 

The origin of the ellipse is given by considering the special 


case in which the particle is going off at m/2 in the barsy, 


25 x = 
; چ٦‎ i 
P 3 cos 6 , = n Es. 


The major axis is the difference between the momentum vectors 
P for the two special Cases with 6" = 0 and 0' < ۲۱ ۱۰6, , 7 p^. The 
minor axis is given simply by 2p'.. 

Computation of the ellipse parameters necessary for construc- 
tion of the ellipse plots was programmed for the IBM 650 by the 
following procedure. 


(a) Calculate the virtual mass of particles No. 3 and No. 5, 
T 7777 د‎ 
A p= 21" › 
ቾ 
where Eis the total center-of-mass energy. 
> a 1 
(b) Calculate p > (max) = [(E, a 2 M] 2, 
where 
fe ች 
E (max) = —, [E “ + MZ - (M 
2E 


2 
፡ 2 + Mg) 1: 


3 
#፡ 
(c) Compute the eilipse parameters for increments of P> in 


the interval O< ٢ É P> (max), using the following expressions: 
1 


ሚጋ : S 3 y 2‏ ل 
کون ۱۳۷۷۷۷۶ E SG‏ 
1 
2 = 
i "ም a - My + M,)] (Bev/c)‏ | ۱ 
ገ E' 20.5 p,[l * (M, - Mg )/ 35]: Bev/c‏ 
where‏ 


z # > ٧ > 
E E UV mas ۷ ۲ ۰و ر‎ 
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In both inelastic cases considered, the negatively charged 
particle isa T , and the positive particle in the final state characterizes 
the event. 

2. Transformation to Center-of-Mass System 

a. Rotation of Chamber Coordinate System 

For the graphical analysis, a transformation of the measured 
momenta and angles of the prongs into the c.m. system is required. 
Propagation of the errors associated with these measured quantities 
into the new reference frame is also necessary. 

Error propagation into the c.m. system was simplified by 
first rotating the measurement coordinate frame so as to bring the 
incident beam direction into coincidence with a transformation axis. 

In the coordinate system of the chamber, System I, Fig. 13, the 
direction vector of an incident T is nearly parallel to the positive y 
axis. System I was rotated to bring the beam direction into coincidence 
with the positive x' axis of System II. 

The beam track in System I has polar and dip angles (¢,; A) 
defined, in terms of direction cosines, by 
Ir 


e 
2 


sinh, = nj), 5 COS) [1 - n 
cos ቀ] = n, /cos زو‎ sind, = n; / sin. 

It can be shown that rigid rotation of the coórdinate System I- 
through an angle of b; about the x axis, followed by a rotation through 
an angle of -A, about the new y axis brings about the desired align- 
ment of the incident track and the +x' axis of System II. Thus, a 
prong 5 represented in System I by 

n-cosAÀcos ií + cos Asin ó j * sin A K 
will, in System II, have the form 

n = sinacos ይ i'* sina 518 ይ |! + cosak', 
with 

sina cos ۵ Cos À cos Ó 

sin a sin P ار ده‎ cos ÀN sin $ 


cos a sinA 
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Fig. 13. Diagram of space coórdinate systems used in error 
propagation to center-of-mass System. (Above) Chamber 


coordinates, System I. (Below) Transformation coordinates, 
System II. 





or 


The real orthogonal matrix M($), Xj) is given by 


COS hy cos ቀ] ር08 ۸۱ sino, sin እ; 
M ($4, X1) - - sind] ር05 ቀ) 0 
- sin À; cos $1 -sin À} sin ڊۉ‎ cos ۸۱ 


In System II the direction of the kth emergent prong is 
specified by 





cosa, = cos A, sin ኢጊ - 81ኬኳእ) cos ኢጊ cos (ቀጊ - ቀ))› 
1 =ተ - 2 
sina, N l -cos a: 
1 1 . 
cos Pi = sina, [cos Xi COS A cos(>, -$))+sin ۸ sin ኢጊ], 
ወ) 
sin By کو‎ a, [cos ኢጊ sin (ቀኒ - ቀ])] ۱ 


It can be shown that the errors in a and B associated with 


errors in $ and À are, for the kth prong, 





2 2 2 2 2 EC 2 2 
0 A oe Are DIE CED 
a lk 2 Pr 2k ۸ 3k እጊ 
2 2 2 2 2 2 2 2 
= + 
VN Bik LZ + (0%, | ር وو م‎ a, 
where the coefficients Ay E. B31 are given by 
Alk = sinB, sin E 
Ak = cos ሰዊ 
Ask = cosh, cos ኢጊ + sinh, sinh, cos($, -$)), 
1 
i” Ge ae 
sin a, 





] | 
E m 35 cos ay sin Py: 


B. = sin à] sin pk 
s sin X 
ዒጊ COSA, 
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b. Center-of-Mass System Quantities 

Following the axis rotation the incident-particle track coincides 
with the + x' axis of System II, and this axis is the transformation axis 
and the reference axis for angle measurements; Or. is the laboratory 
space angle for the kth prong. 

For the kth prong, we can write 


2 


۴ = 


p. cos 0. _ y-n pk 5 ፀጊ 
E n Y E | 
k mn k 


and the center-of-mass momenta necessary for graphicai analysis 


are given by 


* _ 2 EU B دش‎ + 0 
Pk 7 Py [sin 0 * Y (cos 0, - pw l2, (k = 2,3) (V-1) 
k 
where 
p, - "k cos Û, = sina, cos 8 Sino -] ee 
Kun 7> k k k’ kar)? k’ 
Pf FM, 
The angular separation of the two emergent prongs in the c.m. 
> 
system, 6, 4» can be calculated from the invariance of 713 3 › E 
2 D 
DI رد و‎ SD COS, 
MN CENE MR CMM NN 
2 3 ۳۷8" ۰ ۰ 
hence ee 
دګ‎ 1 E, Ez - E, E, + P2 P3 cos 65. 
DOE COS — m * * 9 
۲2 وط‎ 


(V=2) 


with cos 054 = sina, sina, cos (وم - رم)‎ + cos a, cosa, 


and 
* * 2 2 
E. = N Pk + Mi 
From Eq. (V-1) it can be shown that the standard deviation of 


the error associated with the kth prong's momentum (c. m.) is 
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— E 


Pk 2 2 
0 > = — [(C, o ) t (Di. ٣ ) ] , 
Pk p k 
k 
where 
E n M 
k k ر‎ 5 
C, =| — سس‎ y (cos, - 5 ), 
ee De K 7 
/ 
| 2 
Di = Gee [sin 0 cos ج7‎ ) + 85 sin 0], 
and 
1 2r ED 2 27 
E > 
ae 516ا‎ 81 [(sin a, sin BL) ዉቧ + (cos e COS By) os ] ۰ 


The error in 9... a function of the six variables 
(P5: Pz» 4, B5, Q,, ቦን) that define tracks No. 2 and No. 3, was 
obtained by evaluating the required partial derivatives of Eq. (V-2) 
as difference quotients rather than programming in their analytical 
expressions. 

The center-of-mass quantities are, of course, dependent on 
the reaction assumed, i.e., the mass assigned to particle No. 3. The 
computer program first computes the quantities using M, = کک‎ then 
the calculation is repeated with M, = M... 

For each event the two inelastic reactions are considered in 
turn; and the areas, defined by Pa cos 655) and their associated 
errors, are plotted on the appropriate set of curves. For one of the 
two reactions the terminal point of Pa will, ideally, lie (within errors) 
on the ellipse specified by the value of P»: thereby characterizing the 
reaction. A typical event has been plotted in Fig. 1l. The values of 
p, 2 and the standard deviation O x, O x are used in determining 


the best graphical fit. The quantity to be A in each instance 


is 
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One can write EAREN 2 
2 : 3 
1 


2 
(min) ېه‎ o 
The values of کر‎ for the inelastic Interpretations were compared 
for determination of the correct graphical interpretation. For the 
event plotted in Fig. 11, the missing neutral is obviously a neutron; 
ےت‎ = 99. A detailed discussion of the graphical separation of 
inelastic events 1s not given. In summary, separation by this method 
was identical to the separation obtained by the computer program 
discussed in the following section. 
D. Computer Program for Kinematic Analysis 
of Two-Prong Events 

l. General 

Kinematic analysis of two-prong events--i.e., events having 
one incident r track and two charged outgoing tracks--was programmed 
for the IBM-650. The input data consist of the measured values of 
nine dynamical variables x; (specified below) characterizing each 
event, and for each measured variable its associated standard 
deviation, 0. Output data consist of (a) the least-squares adjusted 
value of each of the nine variables required to make a "best fit'' to 
each of three possible interpretations, (e.g., elastic, and inelastic 
with missing = or neutron); (b) the value of se and a number of 
allied statistical consistency tests; (c) the momentum and trajectory 
of the missing neutral particle ms Or neutron). 
2. Choice of Dynamical Variables 

The nine dynamical variables x, (i = 1,2,9) consist of 
three curvature variables (j = 1, 2,3), 

+ 


Kx. ቁ (Bev/c)! 
J Ip, cos Aj 


poe ሦ for incident track, track No. 1 | 
m ; 


-1 for outgoing prongs, tracks Nos. 2 and 3 





Sine 


the three dip tangents 
tan hj, )( 201,2۳ 
and the three azimuthal angles 


ቀ; ۱ 


These variables were chosen because of correspondence with 
quantities most directly measured from the stereo projections on a 
scanning table. Curvature was chosen in preference to the more 
familiar momentum variable to eliminate the following problem: the 
use of momentum as one of the fitting variables introduces an artificial 
mathematical discontinuity at high particle momentum. For example, 
two beam tracks of momentum 4.0 Bev/c and opposite charge are 
indistinguishable by present techniques. Their momenta, however, 
are two large numbers of opposite sign, whereas their curvatures are 
essentially zero in both cases. At lower momenta this difficulty again 
shows up as a skewness in momentum error. 

3. Interpretation Defined 


The general two-prong interaction may be written symbolically 


as 
(M,) +(M,)> (M,) + (M3) + (Mg), 
where 
M, = incident particle (track No. 1) 
M, = struck particle, 
M, = track No. 2 particle, 
M, = track No. 3 particle, 


M, = missing neutral (not used if elastic). 
An interpretation consists, by definition, of two statements: 
1. Whether the event is to be treated as elastic or inelastic; 
2. The set of masses M), M5, M35, M, -- and, for inelastic inter- 
pretations, M,. 


“wma 4፡ሇ”- u 


ክፍ) 





Ec 


4. Kinematic Equations of Constraint 


For any given interpretation the momentum and energy JE, 
(A=1,2,3,4) of the missing neutral, if any, can be written as | 


3 

Fy = (-) 2 cos 6,/K;. (V -3) 
| 3 

F, = (-) 2 sind /K, (V-4) 
3 

F,=(-) 2 يو ورن‎ (V-5) 


Foe eye 2 05533 


4 4 + E] -E,- Ey), (V-6) 


where the energy of the jth ون‎ is ہس‎ from 


(sé A.) > 2 
A > 
J K. J 


with the mass of the jth particle. Here M (Bev) is the mass of 
the struck particle. 

The kinematic constraints can be summarized into two groups 
according to whether the event is treated as elastic or inelastic: 

l. Elastic: Impose four contraints, F, (x) = 0 for A=1,2,3,4. 

2. Inelastic: Impose the single constraint, O(x) = 0 


3 


where Û is defined by the expression 


1 
6 = 1 
M, + lol? Q <0 
where 3 
= D Z 
Q=(-) FE =Fg - a F. 


1 
2 


For Q positive, IQI? is the apparent mass of the missing neutral. 


The quadratic dependence of & on F, makes it possible to 
reach a least-squares fit for which the energy of the missing neutral 


is negative, i.e., F,>0. These erroneous solutions are fortunately 
infrequent and easily detected. 
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The particular kinematic constraints selected above are some- 
what arbitrary. For example, one of the elastic constraints could 
have imposed the coplanarity restriction. Almost any function 
. -Q) such that G(0) was zero could have been used for the inelastic 
constraint. The kinematic constraints cited above were selected first 
for their ease of computation, and second for their essentially linear 
dependence on the dynamical variables in the range of physical 
interest. 

5. Definition of a Kinematic Fit 

In general the measured values xM of the nine variables of 
an event do not satisfy the kinematic constraints. Therefore a least- 
squares adjustment of this input data is required with the measured 
quantities, ca going to the adjusted values, x.: x > x; (121,2-::*9). 

A least-squares fit to a given interpretation is specified by a 
set of dynamical variables 5; satisfying two conditions: 

(a) The constraints are individually equal to zero, i.e., 

F, (x) = 0, for all constraints A = 1,°:: c. For inelastics 
፲ | is set equal to G andc = 1. 
(b) The chi-square quantity, 


9 7 y M È 
70۷ D 
: 68091 ۱ m. 
1-1 1 


is at a minimum (or stationary) point for all variations, 
; Ox., about X; consistent with the kinematic constraints. 
An alternative, but equivalent, formulation using Lagrangian 
multipliers proves more mathematically tractable. This choice 
replaces (a) and (b) by this requirement: 


The modified chi-square quantity 
5 
یں ا‎ 
M(x,,a,)= x +2 E a, F, (x,) (V-7) 


is minimum (or stationary) under arbitrary variations, 6 په‎ Ox. of 


all variables about their stationary values, X.. 


a); 
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6. Statistical Estimates of the Quality of the Fit 
Achievement of a ''fit" 1s a purely mathematical achievement 
ın programming; 1n no way does it assure that the selected interpretation 
is valid. Therefore, a statistical estimate of the quality of the fit 
to a given interpretation is required. 
The adjustment parameters are normalized by introducing the 


quantities 5, 


5, = (x; - x) / ۰. ከክ ን... ۰(۶ 


1 


The foremost quantity for the statistical estimate is, of course, 


2 
M NE 
i-l 
evaluated at the "fit" If x^ is in the doubtful zone, being neither 
unduly large nor sufficiently small, other tests can be quite useful. 
In particular a distinction is made between: (a) "errors"! 

from inherent inaccuracy of measurement and from turbulence, as 
discussed in Section IV; and (b) 'bloopers' caused by miscoding, 
measurement of wrong fiducials, or jumping to a neighboring track, 
or digitizer errors. In addition, it is helpful to distinguish three 
statistical estimates that were used. 
(1) Global: the tests 18.1 < 3. These tests are dependent only upon 
x, at fit, the measured ہو‎ and the standard deviation in the variable, 
9.. 
(ii) Local: the tests | اک‎ ot | < 3. These tests are peculiar to the 
interpretation fitted. Here l g | is the standard deviation in Š; itself. 
Such tests exploit the derivatives of the constraints at fit and are often 
more sensitive than the global test, provided the interpretation is 
Correct. 
(111) Systematic: correlation of curvature and azimuthal angle errors, 
incorrect nominal beam momentum, and neglect of dE/dx corrections 
in curvature measurement. Corrections of this nature were made by 
hand calculation where required for resolution of otherwise ambiguous 


interpretations. 


cd em a r. 
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1. Trajectory of the Missing Neutral Particle 

The variables appropriate to the missing neutral track--1.e., 
curvature, Kg: dip tangent, tan እን and azimuthal angle, $, “are 
calculated for each inelastic interpretation. When a fit has been 


obtalned, ])٥۱:٢٢ 6130107116 5-3 ۲۲۴۰۰ 06 DY 


-1 ۳ " 

K, = (Ps COS ኢጋ) 5 (F, + E, / or (Bev/c) 9 (V -8) 
Lan = F /(F ፡ TE 272 ORE (V -9) 
5 Sw» 1 2 Ss 5" | 
ቀ5 = arctan (F,/F (۰ (V-10) 


where Eas -F,, 
Eqs. (V-3), (V-4), and (V-5). 


-F, is the momentum vector of the missing neutral, 


8. Mathematics of the Fitting Procedure 
a. Linear Constraints 
The procedure introduced here was developed by Frank T. Solmitz 


ent Let us examine the fitting problem for the 


and Horace Taft. 
special case, in which the constraints constitute linear relations 


between the X: 


9 
F(x.) = g, + د‎ F 








x (V-11) 
i=l 
where Fy; = 9F,/9x,. We now require the (9+ c) partial derivatives 
of M, Eq. (V-7) to be zero: 
SM هم‎ E 
m ee 2a, u ااا‎ 
1 ۸ 
OM اس‎ E ۱ ۱ 
s > O: FG.) z 0. (V-13) 


Substituting Eq. (V-12) into (V-13), we get c equations for the a, 5 


M 
= - 14) 
La, Een E (x ^), (V-14 
where a is the symmetric matrix, 
7 
= M ۱ /ሇ..]=5ነ 
Hy (Fy. c (Fj ° د‎ 
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Hence we have 
آتے‎ 


a. = > H F , ነ 1... 
are (V-16) 


which, when substituted back into Eq. (V.12), gives the adjusted values 
OÍ x.. 

i 
b. Nonlinear Constraints 

Solmitz has suggested the adaptation of this method to the 
میں‎ ; . 1 

solution of the fitting problem for nonlinear constraints. : For a point 
x. in the neighborhood of the fitted values X; the nonlinear constraints 


can be approximated by the linear expansion 


۷ V 
Fy (x,) = F, (x; E 0F/0x. | via, -x,; ). (V-17) 


. 1 
1 


The method of successive approximation, using equations paralleling 
Section VD8(a) 1s then used to derive better estimates, x, of the fit 
variables X.. 

The fundamental Eqs. (V-12) and (V-13) carry over directly; 


the quantities Fy. must be here reinterpreted as partial derivatives: 





FS BE A ፡ (V-18) 


Given the vth approximation (x. " 


, a እ) to the solution of the fundamental 
Eq. (V-12) and (V-13), one deduces the next iteration (v+1) froma 


ane V 
first-order expansion in (x; - x, ): 


vtr » v+1 ۷ 
S; E a Fai 0 ور‎ (V-19) 
A 
1 د‎ y (penc 5 sü (V-20) 
EU UU EGER Oc oc ] 

7 ۷ vtl e 
Introducing the two quantities E ‚a, generalizing Expressions 
(V-15) and (V-16), and the modified constraint function 1 , one has 

v _ ነ v v 
EON = 20 a JE ሮ ,), (V-21) 


1 


vtl _ > vl .v 
TT EX I (V -22) 
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vtl 


S. 


; is obtainable by substitution of av" into Eq. (V-19). 


This iteration process can be initiated with ''zeroth-approxi- 


(0) = e l.e., هع‎ - 


mation'' values x, 


9. Statistical Correlation of Adjusted and Measured Variables 


M 
The measured variables x. are assumed to be uncorrelated: 


M. M ` 2 
7 ሦ | = 54 MES , where 7 


0 
is the standard deviation in the measurement of cme However, the 
fitted variables x, are correlated as a result of the fitting process. 
The normalized displacements é (i= 1, ٠٠٠ 9( are also correlated. 

Evaluation of the fit achieved for any interpretation is possible 
only if we establish two references: (a) the standard deviation , 
a of 5. itself; and (b) the average value of x to be expected with 
C constraints, using the fitting procedure. 


These results can be deduced as follows. At x * minimum the 


normalized displacement 5: is given by 


ሠ,” 
and accordingly the variation Ó S; is 
» ET 
02 E Of H , (no): 


| H TEST 
H> V 


Using the result of s lo valid only to the extent that the region 
ےرت‎ : 
near x minimum can be treated as linear, 
Of, Of. SH, 
v pv 


then one deduces 
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Not all standard deviations c ne equal to unity, because 


in a least-squares fitting those dynamical variables least well known 
are, in general, adjusted (shifted) by a greater number of standard 
deviations in the fitting process. The digital computer was pro- 


3 


grammed to evaluate and print out the adjustment variables 5/0 E 


which have in principle a unit- width normal distribution. : 


The mean value of ve can now be established as a special 


result: 


ee 


y? 2% _ dog? 


1 1 


C 


T 2 1 os 
1,۶۶ ١ 


-1 
۷ 


- c, the number of constraints. 


In summary, x should average to | and 4 for the inelastic 
and elastic events, respectively. The 5 distributions for the inelastic 
and elastic events are shown in Fig. 14. 
10. Comparison of Measured and Computed Quantities 

Figure 15 shows a photograph ofa ۲ + ۲ج و‎ + r +n 
interaction. A schematic diagram of the event is shown in Fig. 16. 
This is a typical two-prong event if we consider only tracks Nos. 1, 2, 
and 3. (The kinematic program considers only these tracks.) In this 
event, the missing neutron (neutral track No. 4) strikes a proton, 
track No. 5, and the neutron momentum and trajectory can be calculated 
from the direction and range of this stopping proton. This fortuitous 
circumstance permits the following tests of the kinematic program. 
(a) Does the program, utilizing only measurements of tracks Nos. l, 
2, and 3, achieve a "fit" for the neutron interpretation? (b) Do the 
computed trajectory and momentum of the neutron agree with the 
measured values (based on the vertex, the first bubble of the proton 
track, and the range of the stopping proton)? Table III summarizes 
the comparison between measured and computed quantities. The 


large measurement error associated with the dip angle of tie neutron 
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Fig. 14. (Above) x inelastic interpretation, 250 inelastic 


events 2 "M 
(Below) x elastic interpretation, 250 elastic events. 
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Fig. 15. Photograph of T T+ p > m +m +n event (neutron 
scattering) in the bubble chamber. 
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Fig. l6. Schematic diagram of rT +p +t + r” +n event. 
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is caused by the shortness of the track and the fact that only two points 


were available for measurement. 


x (a typical track) are included in Table IU. 


Table III 


Measured and adjusted quantities 


Value from 
Particle 


error 


0.002 Bev/c 


0.036 Bev/c 
له‎ 
38.0? 


0.025 Bev/c 
(Ro? 
0.5” 


The measurement errors for the 


Measurement Least-squares 


value 


0.267 Bev/c 
93.1; 
و2‎ ۲ 


2299 
2 
0372 6 


74 45 
125092 


Quantity measurements 
proton, range 7۰٦ 
a Ps Greases e 
neutron, Pa 0.256 Bev/c 
No. 4 ቀ; 95.19 
0 
Ay 19 ت3‎ 
length 05825 cm 
(meas.) 
5 2 
u X el 
2 
Xn 
Tt Pz 0.388 Bev/c 
No. 3 $4 74,19 
O 
Az 11.9 
length 13.138 cm. 


(meas. ) 
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VI. CLASSIFICATION OF EVENTS 
All events were first examined for the value of the coplanarity 


function ۴ and catalogued as coplanar Ina < 4) or non- 


coplanar (/o 24). 

The noncoplanar events were considered to be definitely in- 
elastic, and were subdivided as follows: 

(A.1) Strongly separable by کی‎ i.e., r نات‎ x“ > Ser 
vice versa. Final classification of these events corresponds 
to the smaller "e 

(A.2)  Tentatively ambiguous, i.e., both x^ 0 and x^. <5. These 
events were examined for delta rays and ionization of the 
positive track, and kinematically examined on the Po? Py 
plot described below. 

(A.3) All others - not having satisfactory for either the ፐ" or 
neutron interpretation. These events were examined on the 
scanning desk for measurement and (or) coding errors, 
remeasured or adjusted by template fittings, and finally 
forced into category A.l or A.2. 

The coplanar events, لكب ا‎ < 4, were divided into classes 
in an analogous way: 

(B.1) Strongly separable as elastic, by X 

(B.Z) Elastic or Tr (ambiguous); neither ionization nor delta rays 
can possibly resolve the ambiguity. 

(B.3) Elastic or neutron (ambiguous); possible to resolve by 
identification of positive track, on the basis of ionization and 
delta rays. 

(B.4) Neutron or nO (ambiguous); separable as class A.2. 

(B.5) No fit achieved; events remeasured. 

All the events with | E /o اب‎ < 4 and fitting only inelastic 
interpretations were re-examined for gross measurement errors. 
These errors can induce a spurious mÜ or neutron fit. 

Positive-particle identification, based on ionization, was used 
in the classification of otherwise ambiguous events. The incident m 
was used as the local minimum-ionization standard. (Calibration of the 


chamber for bubble density as a function of ionization was not feasible 
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for two reasons: (a) chamber sensitivity was kept high to detect neutral- 
particle end points; (b) superposition of bubbles--especially for steeply 
dipping tracks--makes bubble counting unreliable.) Low-momentum 
protons were readily identified in the chamber. 

It is quite possible for an inelastic event to be kinematically 
ambiguous, even in the limit of exact momentum and angle measurements. 
To see this let p, be the observed momentum of the positive prong, and 
let Po be the calculated momentum of the missing neutral (e or neutron). 
Normally the energy available to the positive and neutral particles, as 
calculated from the conservation of energy and the known masses of 
the incident and struck particles, allows assignment of only one mass 
to the positive particle and one mass to the neutral particle (either 
proton mass and nm mass, or m" mass and neutron mass). However, 
to the approxımation that the small rs T, n-p mass differences can 
be neglected, kinematic separation fails for Po = P,- Although such 
events are kinematically ambiguous, they may, of course, be classified 
on the basis of delta rays or ionization along the positive track. 

Figure 17 shows the pg, p, distribution of a group of 26 events 
of class A-2. The ambiguous Po = Py line is broadened by the typical 
measurement error of pj: the error increases with increasing positive 
momentum. For the most part, those events with good 20. separation 


have the better interpretation (m or n) outside the ambiguity region. 
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Fig. 17. Scatter diagram of positive and neutral particle 


momenta (26 events). Kinematically ambiguous region 
delineated. 
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Vil. RESULTS 
A. Cross Sections 

The total number of interactions, corrected for over-all 
scanning efficiency and detection dependence upon the orientation of 
the plane of emergent prongs is 1430 + 43. After correcting for the 
small-angle elastic scattering, the total number of events occurring 
in the acceptance volume of the chamber is increased to 1452 + 48. 

The total number of interactions and the integrated pion path 
length in the acceptance volume provides a direct measurement of 
the total cross section. Beam tracks crossing the acceptance volume, 
as projected in one view, were counted in every tenth frame of film. 
Based on this beam count and the average track length across the 
acceptance volume, a total pion path length of ” cm+ 1.1% 
was computed. This length includes corrections for muon contam- 
ination (1.0 + 1.0%) and electron contamination (1.8 + 1.2%). i: The 
T -p total cross section found is 53.3 + 2.4 millibarns. The cross 
sections for all zero-, two-, four-, and six-prong and V-particle 
events observed are listed in Table 1V, together with their relative 
frequencies. 

Table V gives the relative frequencies and cross sections for 
the elastic-scattering and single-pion-production reactions. 

B. Elastic Differential Cross Section 

Figure 18 shows the elastic T -p differential cross section. 
This differential cross section displays the narrow diffraction peak 
in the forward direction, and after dropping to a distinct mimimum 
from 66? to 108? (center-of-mass system) shows a pronounced hump 
in the backward hemisphere. The backward hemisphere contains 
35.6% of the events. The differential cross section at 0°, do (0)/6Q, 
extrapolated from the curve is 15 + 4 mb/sterad. Cool, Piccioni, 
and Ci give a value of 9.4 (center-of-mass system) on the basis 
of a 46-mb total cross section. Adjustment of the imaginary part of 
the forward-scattering amplitude for the larger total cross section 
of 53.3 mb yields a prediction of 12.6 mb for the elastic differential 


0 
cross section at Ü . 
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Table IV 


Phenomenological classification 


Fraction of Cross section 
Reaction No. of events total events (mb) 
Zero prong ئ2‎ 9 0.176 + 0.006 9.4+0.3 
Two-prong | E aa 0.814 + 0.024 43.4 + 1.3 
Four-prong 1222 0.008 + 0.001 0.4 + 0.001 
Six -pr ong 0 0 0 
Associated 3-21 0.002 205001 O.1 = 0-001 


production 


$ * 
Includes 22 + 5 events in the interval 1 >cos o - 20.95 ٩ 6 


differential cross section, Fig. 18. 


Table V 


Elastic and single-pion production 


No. of events Fraction of Cross section 
Reaction analyzed analyzed events (mb) 
Elastic 261 0.483 21.0 + 1.1 
ና کن‎ 80 0.148 6.4 
ሸ ተ ተክ 199 0.369 16.0 
All other — — 90 


interactions 


Total cross section = د د5‎ 3 2 4 
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Fig. 18. Elastic m-p differential cross section at 0.939 Bev/c. 
Events corresponding to scattering of less than 189 in the 
center-of-mass system are omitted. 
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C. Distribution of Charged Particles in Single-Pion Production 

The momentum distributions of the "T and proton from 
T +p>m +p+ ne reactions are given in Fig. 19. The momentum 
distributions are compared with those predicted by phase-space 
considerations. The average momentum of the observed m is some- 
what higher than the average value predicted by phase-space arguments, 
and are consistent with the momentum distribution anticipated from the 
isobar model. 0 

In Fig. 20, the momentum histograms of the charged pions 
1۱ ۳ the reactionT + ۳ —>wT + l ፒ0 Again the momentum 
distributions are compared with the phase-space distributions. A 
marked departure well outside statistics is seen; the o peaks at a 
momentum of 125 Mev/c greater than center momentum of the 7 
distribution. 

No attempt was made to identify multiple-pion-production events 
during this experiment. The number of four-prong events observed 
(0.8 + 0.1% of total events) indicates that pion multiplicities were 
infrequent at this value of incident n momentum. The single-pion- 
production events analyzed in this experiment provide a reference 
point for analysis of single-and multiple-pion production at higher 
values of momenta, specifically the 1.12 and 1.24 Bev/c data now in 


the initial stages of analysis. 
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Fig. 19. The momentum distributions of the protons and T's 
trom the m pn a mÜ reaction. Expected phase- 
space distributions are shown for comparison. 
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Fig. 20. The momentum distribution of m's from the 


T^ t p — T +T" +n reaction. Phase-space distributions 
are shown. 
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